Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 25 February 2010 



(MN style file v2.2) 



Jeans Instability of Palomar 5's Tidal Tail 



o 

(N 

in 

(N 

u 

o 



> 

o 

00 
O 

o 



Alice C. Quillen & Justin Comparetta 

Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA; 
aquilleniSpas.rochester. edu; jcompare@pas. rochester. edu 



25 February 2010 



ABSTRACT 



Tidal tails composed of stars should be unstable to the Jeans instability and 
this can cause them to look like beads on a string. The Jeans wavelength and tail 
diameter determine the wavelength and growth rate of the fastest growing unstable 
mode. Consequently the distance along the tail to the first clump and spacing between 
clumps can be used to estimate the mass density in the tail and its longitudinal velocity 
dispersion. Clumps in the tidal tails of the globular cluster Palomar 5 could be due 
to Jeans instability. We find that their spacing is consistent with the fastest growing 
mode if the velocity dispersion in the tail is similar to that in the cluster itself. While 
all tidal tails should exhibit gravitational instability, we find that clusters or galaxies 
with low concentration parameters are most likely to exhibit short wavelength rapidly 
growing Jeans modes in their tidal tails. 



1 INTRODUCTION 



A thin self-gravitating cylind er of gas can clump into a series 



of bead like structures (e.g., Chandrasekhar &: Fermil 1 19531 : 



IOstrikei|[l963 : lElmegreenHl979l . Il994l ') . This model has been 



use to propose explanations for clumps of star formation 
and s tar clusters along spiral arms (E lmegreen & Elmcgrccn 



1983 ) and in tidal tai l s of galaxies (lElmegreen fc Efremovl 
1999 : IDuc et al.l I2OO0I : ISmith et all I2OI0I '). The sausage or 
varicose instability is a rough analogy to the Plateau- 
Rayleigh instability where surface tension causes a stream 
of water to become a series of droplets. A cylinder of 
stars can also contract into a series of clumps. In this 
case the stability has been c alled the Jeans instability 
iFridman fc Polvachenkol 1 19841 ) even though this term is 
more often used to describe the growth and decay of plane 
waves through a homogeneous stellar medium. The instabil- 
ity for a rotating stellar cylin der has been investigated by 
iFridman fc Polvachenkol (| 19841 ') (see their Chapter 2). In the 
limit of infinitely slow rotation their analysis can be used for 
the non-rotating cylinder. 

Nearly equally spaced clumps observed i n the tidal de- 
bris of Palomar 5 l|Odenkirchen et afllioOll . 120021 ') has not 
been interpreted in terms of the Jean instabili ty but instead 
prim arily in terms of oscillations in the cluster (|Gnedin et al.l 
I1999I ') ca used by a previous passage through the G alac- 
tic disk (|Odenkirchen et al.|[200^ : iDehnen et al.ll2004h . Al- 
ternative explanations for structure in tidal tails include 
the effec t of structur e in the dark matter halo, as ex- 
plored bv lMaver et al.l (^2002) and variations in density along 
the tail caused by epicyclic motion s in the stellar orb its 
l|Kiipper et al.ll2008l : [just et al.|[2009l : iKiipper et al.ll201(]| ). 

Here we consider the possibility that clumps in a tidal 
streams could be caused by gravitational instability. 



We review the properties of the Palomar 5 clust er and 
tail s ystem. Its distance is estimated at 23.2 kpc (|Harrisl 
Il996l l. The ve locity dispersion in the ta il is small, measured 
at 2-4 km/s (|Odenkirchen et al.l |200S^ ). The dispersion in 
the cluster itself is small, ^ 0.3 km/s (correcting for bina- 
ries |Odenkirche^^eF^l][2QQ2)5 s^nd the cluster is suspected to 
have low concentration parameter, Ct, in the range 0.4-0.8 
lOdenkirchen et al.|[2003 : lDehnen et al.ll2004l '). The width of 
the tail is 120 pc and the tidal radius of the cluster half of 
this (Odcnkirchcn ct al. 2002, 200j). The rn ass of the cluster 
is about 5000 Mq (jOdenkirchen et al.ll200^ ') and that in the 
tail exten ding 10° exceeds the mass in th e cluster by a factor 
of a few (|Odenkirchen et al.|[20o3 . |2003|). The tail has been 
detected up to 22° from the cluster (|Grillmair fc Dionatod 
I2OO61). The linear mass density in the two tails as a function 
of distance along the ta il drops slowly as a funct ion of dis- 
tance from the cluster (|Odenkirchen et al.ll2003|) as would 
be expected from a tida l evolution model ( Johnston et al.l 
1 19991 : iDehnen et aT]|2004l ') . 



Clumps are most easily seen on the southern side of 
the tidal tail with distance to the first clump about 1.7 de- 
grees from the cluster center and that between the first 
and second clump abo ut 2 degrees (using Figure 14 by 
lOdenkirchen et al.|[20o3 that also shows that extinction can- 
not account for the density variations). At a distance of 23.2 
kpc, 2 degrees corresponds to 800 pc and we can treat this 
as the wavelength of a possible growing unstable mode. It 
is useful to take the ratio of the wavenumber to the tidal 
tail radius, ro, (diameter/2). We estimate kro ~ 0.5 using 
To ~ 60 pc. 
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2 JEANS INSTABILITY OF A STELLAR 
CYLINDER 

We briefly review the gravitational instability of a stellar 
cylinder to varicose or sausage like compressive modes of os- 
cil lation. The dispersion relation i s derived more rigorously 
bv lFridman fc Polvachenkol (| 19841 ') (in their Chapter 2) and 
covering the case of a homogenous rotating stellar cylinder. 
We can adopt the simplest assumption of a constant linear 
mass density, ^o, but a deformed boundary of radius R{z) 
and a corresponding interior density, p, that is independent 
of radius interior to R. We are working in cylindrical coor- 
dinates with z oriented along the cylinder. To first order in 
a perturbation amplitude a 



p{z) = po (1 + d cos fcz) 



R{z) = ro 



^1 — — cos kz^ 



(1) 



Here we have assumed that our perturbation amplitude, a, 
is small. The relation between the perturbation amplitude 
in radius R and density p is determined to first order by the 
assumption that the linear mass density does not depend on 
z and /io = PottTq. 

To estimate the perturbation to the gravitational po- 
tential we can consider the potential perturbation caused 
by a massive wire with mass density p(r, z) = S{r = 
0)(/io -I- pi cos(fcz)). The part of the gravitational potential 
fluctuating with z at a distance r from the wire and at z is 



$(r, z) 



Gpi 



cos{kz') 



+ (z - z')^ 



zdz' 



-2Gpi cos{kz)Ko{kr) 



(2) 



where Ko is a modified Bessel function of the second kind. 

To estimate the gravitational potential perturbation we 
must integrate the Bessel function over radius interior to 
the cylinder boundary. For kr < 1 the asymptotic limit 
Ko{kr) — > — ln(fcr/2) — 0.5772.. where the constant is the 
Euler-Mascheroni constant. Integrating equation [2] out to 
our boundary R (given in equation [1} the potential pertur- 
bation 



<E>i(2,r- = 0) 



AirGpoa cos kz 



In 



0.5572 



rdr(3) 



Deviations on the boundary only contribute to second order 



"l>i(z) « — 2G/ioaj 1 — In kro \ cos kz 



(4) 



where we have used a constant of 1 inside the absolute value 
as it is a more accurate match to the integral of the Bessel 
function than the constant 0.62 given by the asymptotic 
limit. To first order in a the same relation would have been 
estimated if we had assumed no boundary deformation but 
a longitudinal variation in linear density p. This situation 
corresponds to longi tudinal modes rather than sausage or 
varicose modes (e.g.. IOstrikei|[l964l ). 

We now consider density and potential perturbations 
that are proportional to e'('=^-'^*>. We assume a phase space 
stellar distribution function in the cylinder 

/(2,«,^) = po[l + ac'(-*-'=^']^(^;) (5) 

consistent with our density perturbation (equation [T]) . Here 
fa{v) is a Gaussian velocity distribution and we need only 



consider the distribution as a function of the z velocity com- 
ponent. 



fa 



: exp 



2^ 



(6) 



where ct* is the stellar velocity dispersion in the z direction 
We use a linearized version of the coUisionless Boltz- 
mann equation, 

d^i dfo 
dz dv 



d d 
di+^d-z 



fi 



= 



(7) 



where we have denoted the unperturbed and first order per- 
turbation to the distribution function as /o and /i (see equa- 
tion [5}. Using equation |4] for the potential perturbation we 
find 



2Gpo\l-Hkro)\ 



dv- 



kv 



exp 



^2-Kal ^-kv 
valid for fcro < 1. We can rewrite this as 

udu 



1 



In fcro I 



1 e" 







where Z — a;/(\/2A:a,) and 

= - 2 

^ ~ 2Gpo " (fcjro)2 



(8) 



(9) 



(10) 



where k. 



J 



''-^^Po jg ^]jg jgans wavenumber and kj = 



2n/\j with A J the Jeans wavelength. 

The value q is analogous to the Toomre Q instability 
parameter. The dispersion relation is similar to that for 
the Jeans instability for a plane wave propagating through 
a homogenous stellar medium. The above dispersion re- 
lation that we have estimated roughly from the poten- 
tial perturbation at the center of the cylinder is consis- 
tent with that in equation (14 ) more rigorously derived by 
Fri dman fc Polvachenkol (|l984[ ) for the homogenous rotating 
cylinder. 

Assume that u — i"/ leading to unstable solutions. Then 
we can integrate equation [9] finding 



In kro I 



(erf7' - l) + 1 







(11) 



where 7' = 7/(\/2fca,). This dispersion relation implies 
that a "long cylinder is unstable to any large but fi- 
nite thermal scatter of particl es in longitudinal velocities" 
(jPridman fc Polvachenkall984j just before their equation 21; 
also see their Figure 10 showing a numerically measured sta- 
bility boundary for the non-rotating cylinder with numerical 
calculations done by S. M. Churilov). 

The above equation gives a relation between growth rate 
and wavelength. As 7' depends on k the growth rate is most 
easily discussed in terms of 7'fero = 'yra/at or the growth 
rate in units of the tidal tail's radius and velocity dispersion. 
Interestingly the growth rate seems to be of order a,/ro 
rather than of order y^Gpok. However if g ~ 1 and kvg ~ 1 
then the two growth rates are similar. 

In Figure [1] we show growth rate versus wavenumber for 
different values of q. We see that at each q there is a mode 
that is maximally unstable or grows fastest. In this flgure 
we have plotted 7ro/cr* versus kro so that the growth rate 
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Figure 1. Growth rate vs wavelength for 5 different values of 
q. The wavelength of maximum growth rate depends on g with 
smaller wavelength (larger wavenumber k) having a faster growth 
rate if q is smaller. The x axis is the wavenumber times the tidal 
tail radius. The y axis is the growth rate in units of ro/cr*. For 
Palomar 5 we have estimated from the distance between clumps 
that the fastest growing mode has krg ~ 0.5 corresponding to 
q ^ 0.9 and a growth rate of fro/at ~ 0.2. The estimated growth 
rate and wavenumber of the fastest growing mode for Palomar 
5's tails is shown with a solid circle in the middle of the Figure. 



(shown on the y axis) is independent of the wavenumber. 
At larger k, the growth rate drops to zero and there is no 
longer instability. However for any q, even those greater than 
1, there are wavenumbers that are unstable. As q increases 
these wavenumbers correspond to increasingly long wave- 
lengths and the growth rate becomes slower and slower. For 
g < 1, larger k wavenumber (or shorter wavelength) modes 
can grow. 

For sufficiently small q we must keep in mind that we es- 
timated the gravitational potential for kro < 1. For kro > 1 
(corresponding to wavelengths shorter than the cylinder di- 
ameter) the potential perturbation should be proportional 
to k~^ and we should recover the Jeans instability for plane 
waves. The condition g < 1 is equivalent to the Jeans wave- 
length being shorter than about twice the diameter of the 
stream (see equation 1 10|) . For the Jeans instability of plane 
waves, the larger the wavelength, the faster the growth rate. 
Consequently when q <^ 1 we expect that the fastest grow- 
ing wavelength will be of order the one that fits within the 
diameter of the tail, at which point the dispersion relation 
deviates from that estimated for a plane wave. 



3 APPLICATION TO PALOMAR 5 

In the introduction we estimated from the distance between 
the clumps on the southern tail of Palomar 5 that the fastest 
growing mode has kro ~ 0.5. Using Figure 1 we estimate 
that q ~ 0.9 in Palomar 5's tail as that value gives the cor- 
rect wavelength for the fastest growing mode. As q depends 
on both linear mass density and longitudinal velocity dis- 
persion we have a relation between these two quantities. We 
first review estimates of the linear mass density based on 



observations and then we place constraints on both these 
quantities. 

The tidal tail linear mass density for a cluster that is 
losing mass at a rate M is 



dM 
dz 



M 



dt 
dz 



(12) 



-Mrt 



where ^ is the rate that stars are drifting away from the 
cluster and z gives distances along the tail. 
At a radius R from the Galaxy center, 

dz dO. 

Itt ~ Ir^ 

where S1(_R) is the angular rotation rate of a particle in a 
circular orbit in the galaxy radius R. This is approximately 
consistent with the m ore precise semi-analytical model by 
Ijohnston et al.l (jl999l ). For a galaxy with a fiat rotation 



(13) 



dz 



(14) 



For a tidal radius of 60 pc, a galactic rotational velocity of 
Vc = 220 km/s and distance of 23.2 kpc, the drift rate in 
Palomar 5's tidal tail is dz/dt ~ 0.6pc Myr~^. The above 
estimate is only approximate as a better estimate would 
take into account the cluster orbit and adopt a more re- 
ali stic Galactic potential (e .g., as exp lored in the appendix 
by lOd enkirch en et al.ll2003l: a lso se e Ijohnston et al. I ll999l : 
lJuste t al. 20091: iKiipper et aLllioioh . 

Using the estimated drift rate of 0.6 pc Myr~^ 
the mass loss rat e of M ~ 5Mq Myr~^ estimated by 
lOdenkirchen et al.l (|2009l ) corresponds to a linear density of 
H ~ 8Mq pc~^. This is three times lower than a linear mass 
density estimated from stellar number counts, with the num- 
ber counts in the tails over 10 degrees exceeding that in the 
cluster by a factor of ab out 2 and the estimated total clus- 
ter mass of ~ 5OOOM0 (|Odenkirchen et al.|[20o3 '). We note 
that mass segregation may affect estimates of the tail mass 
density based on star counts (,Koch et al.,i2004 ). 

We can write for q 



2Gu ^ ^ ( 



0.3 km s 



^J■o 



IOM0 pc-i 



(15) 



where we h ave inserted the estimate d velocity dispersion of 
the cluster (jOdenkirchen et al.ll2002l ) and a scaling value for 
the linear density. We can see that q may be of order 1 for 
Palomar 5's tail just by looking at the above scaling values. 

The position of the first clump gives us an estimate of 
the growth timescale for the instability. The first clump is 
visible only about 2 degrees away from the cluster center. 
This corresponds to 800 pc or a time 1.3 Gyr using the 
above drift velocity. This implies that the growth rate is 
approximately 7 ~ 0.7 Gyr~^ (assuming a single exponential 
growth timescale to the first clump). 

The growth rate is related to the longitudinal velocity 
dispersion in the tail so the estimated growth rate gives an 
estimate for the stellar velocity dispersion in the tail. For 
q ~ 0.9 we expect that the growth rate has 7^ — 0.2 (using 
Figure 1) consequently we can estimate 



a, ~ ~ 0.2km s 



(16) 



This dispersion and our previous estimate for q implies that 



4 



the mass density in tiie tail is about fj, ~ 4Mqpc~^ . This is 
in between the density estimated from the mass loss rate, 
/X ~ 8M0pc~^, and that estimated from star counts ~ 2 — 
3M0pc"^ Had we used a lower value of the constant inside 
the absolute value in equation |4] we would have estimate a 
slightly higher growth grate and linear density. 

The velocity dispersion we estimate above is sig- 
nificantly lower t han the 2-4 km/s measured by 
lOdenkirchen et alj l|2009l ). We note that if the dispersion 
is isotropic a velocity dispersion of 1 km/s would cause the 
tidal stream to spread at a rate of about 1 kpc Gyr~^. This 
would exceed the width of the observed tail a few degrees 
from the cluster (taking into account the estimated drift 
rate). Consequently a velocity dispersion in the tail is ex- 
pected to b e similar to that estima ted for the cluster or 
~ 0.3 km/s (jOdenkirchen et al.l[20oi ) and so consistent with 
our estimate. Numerical simulations bv lDeh ncn ct al. (2004) 
have also concluded that the tail is "cold." iKiipper et al.l 
suggest that stars escape the cluster from Lagrange 
points and so the tail could have velocity dispersion lower 
than the cluster dispersion. In short we find that the tail 
density and velocity dispersion could be consistent with 
an interpretation of clumps in Palomar 5 as the fastest 
growing Jeans unstable mode of a stellar cylinder. Here 
we used a distance of 2 degrees between clumps to esti- 
mate the wavelength of the fastest growing mode but there 
could also be structure o n 1 degree scales (see Figure 14 by 
lOdenkirchen et al.ll20o3 ). We used a 2 degree scale because 
our dispersion relation is appropriate for fcro > 1. Future 
work could improve calculation of the dispersion relation 
for the fcro ~ 1 regime. 



4 ESTIMATING THE TAIL DENSITY 
DURING TIDAL STRIPPING 

The mass loss rate due to tidal disruption depends on the 
radial density distribution of the cluster and on the Galactic 
tidal field. Thus tidal evolution models should be able to pre- 
dict the linear mass density and velocity dispersion in tidal 
tails. To estimate these quantities we explore the simplest 
model, that where both background galaxy and cluster are 
described as spherically symmetric isothermal spheres with 
dispersions Ug and ad, respectively. It may also be useful to 
write things in terms of velocities of particles in circular or- 
bits Vg — \/2(Tg . A balance of tidal force versus self gravity 
gives a relation for the tidal radius rt of the cluster 



dM 
dz 



= = 
dz 



G 



(19) 



R 



(17) 



where the cluster is at a radial distance R from the galaxy 
center. This gives a mass loss rate for the cluster 



(18) 



Material stripped from the cluster at the tidal radius will 
be on orbits that move away from the cluster because they 
have a different angular momentum. We let R — bvg with 
unit-less parameter b that depends on the cluster orbit with 
6 = if on a circular orbit and b of order unity for a radial 
orbit. Using equation (|12|l and (|14|) and the above relation 
for Af , we can estimate the mass density in the tidal tails as 



The dispersion in the tail should be similar to the veloc- 
ity dispersion of the cluster, a^, and so using our estimate 
for /I, 



_l_ 

4b' 



(20) 



The above value for q is of order unity for a moderately 
eccentricity orbit. We expect that streams from tidally dis- 
rupting isothermal spheres with a moderately eccentric orbit 
have q of order 1 and so are likely to have fastest growing 
mode with wavelength of order a few times the radius of 
the stream. In this regime the growth rate can either be 
estimated with a velocity ^/Gjl or a velocity ad- 

The estimates in this section are for an isothermal 
sphere which has radial density profile p cc r~^. This can 
be compared to the core of a King model that has constant 
density. The concentration of the cluster can be described by 
a concentration parameter Ct = \ogj^Q{rt / Ro) where Ro is the 
core radius and rt the tidal radius. This parameter is used 
to describe King models. Exterior to -Ro, clusters with high 
concentration parameter have density dropping more steeply 
with radius than those with low concentration parameter. If 
the density drops less steeply with radius than an isothermal 
sphere, as would be true for a cluster with low concentra- 
tion parameter, then we would expect that the linear mass 
density in the tidal tail is higher than predicted from an 
isothermal sphere. This in turn we expect would lower the 
q of the tail, allowing shorter wavelength modes to grow. If 
the cluster or galaxy is centrally concentrated then q would 
be high, only very long wavelength modes would grow, and 
they would grow very slowly. 

In summary we find that tidal streams from disrupt- 
ing isothermal spheres have q of order 1, with lower and 
higher concentration disrupting objects having tails with 
higher values and lower values of q, respectively. We ex- 
pect the separation between clumps to be set by the fastest 
growing wavelength which we expect should be of order a 
few times the tidal radius of the tidally disrupting cluster. 
If the cluster is small then the timescale for stars to move 
away from the cluster or galaxy is long. But in this case the 
instability takes longer to grow because the density is lower. 
The two effects cancel out to leave an e-folding length scale 
for the growth of the instability in the tidal-tail that is only 
dependent on the tidal radius. 



5 SUMMARY AND CONCLUSION 

We have investigated Jeans instability of a stellar cylinder 
in th e context of tidal tails. As had previous works found 
(e.g., iFridman fc Polvachenkcl Il984l ). we find there are al- 
ways modes that can grow. The wavelength and growth 
rate of the fastest growing mode depends on a g parame- 
ter q — 2G'fj,o ~ 2/(fcjf'o)^ that is analogous to the Toomre 
Q parameter. Here p is the linear mass density, a, is the 
longitudinal velocity dispersion, ro is the tail radius and 
fcj the Jeans wave- number. For g ~ 1, the fastest growing 
mode has wavelength a few times the diameter of the tail 
and growth rate approximately equal to the tail velocity 
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dispersion divided by the tail diameter. Because the wave- 
length and timescale of maximum growth rate depend on 
the q value and velocity dispersion in the tail, both can be 
roughly estimated from the spacing of observed clumps in 
the stream and the location where they appear. 

For a cluster or galaxy that is approximately an isother- 
mal sphere we find that the disrupting system will have tidal 
tails with q of order 1. Clusters or galaxies with low con- 
centration parameters and on eccentric orbits could have 
tidal tails with lower q parameters and so their modes of 
maximum growth rate should have shorter wavelengths and 
faster growth rates. More concentrated objects would have 
tails for which the mode of maximum growth rate is more 
slowly growing and longer wavelength than for low density 
or low concentration parameter objects. All tidal tails could 
exhibit Jean instability, however if the cluster or galaxy is 
sufficiently concentrated then the wavelength of the fastest 
growing mode in the tail would be large and growth rate suf- 
ficiently long that the instability may be difficult to detect. 
The opposite is true for objects that are close to being com- 
pletely disrupted; Jeans instability should be particularly 
prominent. 

The tidal tails of Palomar 5 exhibit clumps with spacing 
approximately a few times the tail width. From this we have 
estimated a tail density and find it approximately consistent 
with observational estimates. This suggests that the clumps 
in Palomar 5's tails are consistent with the fastest growing 
Jeans unstable mode. Our estimated veloc i ty dis persion is 
lower than measured bv lOdenkirchen et al.l l|2009l ) but con- 
sistent with the estimated cluster velocity dispersion. Palo- 
mar 5 has a low concentration parameter and so its tails 
would be expected to exhibit a quickly growing, short wave- 
length, fastest growing Jeans unstable mode. 

We find that tidal tails from disrupting clusters or dwarf 
galaxies should display gravitational instabilities and from 
measurements of the wavelengths of the fastest growing 
mode, properties of the tail can be directly measured in a 
fashion only weakly dependent on the orbit or host galaxy 
dark matter profile. More detailed understanding of the in- 
stability may yield constraints on the orbit and tidal mass 
loss rate that are complimentary to those from studies of 
shape of the tail on the sky or its velocity gradient. 

Here we have ignored the nature of the stellar orbits 
in the tails as they move through the galaxy. Ifowever den- 
sity variations will be induced in the tail by its orbit. For 
Palomar 5 the growth timescale for the instability exceeds 
the orbital period so it may be sufficient to consider the tail 
density averaged over the orbit. In other systems the orbit 
period may be shorter than the growth timescale for the 
fastest growing mode. In this case at orbit apocenters a lin- 
ear density increase could cause shorter wavelength modes 
to grow. In both cases as the tail expands and density drops, 
shorter wavelength unstable modes may become stable. 

The Jeans instability m odel for clumps in t idal tails 
differs from that proposed bv lKiipper et al.1 l|2008h who in- 
terpreted clumps in the tidal tails of Palomar 5 in terms 
of epicyclic over-densities. By studying systems of different 
densities and with different orbits, future studies should be 
able to differentiate between the possible explanations for 
clumping in stellar tidal tails. Future work will determine 
if the clumping seen in N-body simulations of higher den- 



sity tidal tails from dwarf galaxies jComprretta fc Quillen I 
is due to Jeans instability. 
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